In [2]:
%matplotlib notebook
from charistools.modelEnv import ModelEnv
from charistools.hypsometry import Hypsometry
import matplotlib.pyplot as plt
#pylab.rcParams['figure.figsize'] = (16.0, 8.0)
/Users/brodzik/miniconda/envs/charistools/lib/python2.7/site-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.
warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.')
In [3]:
tileConfigFile = '/Users/brodzik/charis_training_2015/' + \
'module05_temperature_index_model/modis_tiles_config.ini'
topDir = '/Users/brodzik/projects/CHARIS'
myEnv = ModelEnv(topDir=topDir, tileConfigFile=tileConfigFile, verbose=True)
charistools.modelEnv: read MODIS tile configuration from /Users/brodzik/charis_training_2015/module05_temperature_index_model/modis_tiles_config.ini
charistools.modelEnv: model_top_dir changed to /Users/brodzik/projects/CHARIS
charistools.modelEnv: read MODIS tile configuration from /Users/brodzik/charis_training_2015/module05_temperature_index_model/modis_tiles_config.ini
In [12]:
import os
os.getcwd()
Out[12]:
'/Users/brodzik/charis_training_2015/module05_temperature_index_model'
In [3]:
station_tempFile = '/Users/brodzik/projects/CHARIS/forcing_data/Observed/air_temperature/derived/' + \
'v1.0/hunza_observation_tavg_profile_100m.2001.dscale.v1.0.asc'
drainageID = "IN_Hunza_GDBD"
year = 2001
uncorrected_tempFile = myEnv.hypsometry_filename(drainageID=drainageID, type='temperature_by_elevation', year=year)
corrected_tempFile = '/Users/brodzik/projects/CHARIS/forcing_data/Downscaled/ERA-Interim/day/' + \
'corrected/by_elevation/IN_Hunza_GDBD.2001.0100m.corrected.v2.temperature_by_elev.txt'
print station_tempFile
print uncorrected_tempFile
print corrected_tempFile
/Users/brodzik/projects/CHARIS/forcing_data/Observed/air_temperature/derived/v1.0/hunza_observation_tavg_profile_100m.2001.dscale.v1.0.asc
/Users/brodzik/projects/CHARIS/forcing_data/Downscaled/ERA-Interim/day/uncorrected/by_elevation/IN_Hunza_GDBD.2001.0100m.v1.temperature_by_elev.txt
/Users/brodzik/projects/CHARIS/forcing_data/Downscaled/ERA-Interim/day/corrected/by_elevation/IN_Hunza_GDBD.2001.0100m.corrected.v2.temperature_by_elev.txt
In [4]:
station_temp_hyps = Hypsometry(filename=station_tempFile)
uncorrected_temp_hyps = Hypsometry(filename=uncorrected_tempFile)
corrected_temp_hyps = Hypsometry(filename=corrected_tempFile)
In [ ]:
#%ls /Users/brodzik/projects/CHARIS/forcing_data/Downscaled/ERA-Interim/day/corrected/by_elevation
In [5]:
print(station_temp_hyps.data.columns)
print(uncorrected_temp_hyps.data.columns)
print(corrected_temp_hyps.data.columns)
Index([u'1450', u'1550', u'1650', u'1750', u'1850', u'1950', u'2050', u'2150',
u'2250', u'2350', u'2450', u'2550', u'2650', u'2750', u'2850', u'2950',
u'3050', u'3150', u'3250', u'3350', u'3450', u'3550', u'3650', u'3750',
u'3850', u'3950', u'4050', u'4150', u'4250', u'4350', u'4450', u'4550',
u'4650', u'4750', u'4850', u'4950', u'5050', u'5150', u'5250', u'5350',
u'5450', u'5550', u'5650', u'5750', u'5850', u'5950', u'6050', u'6150',
u'6250', u'6350', u'6450', u'6550', u'6650', u'6750', u'6850', u'6950',
u'7050', u'7150', u'7250', u'7350', u'7450', u'7550', u'7650', u'7750'],
dtype='object')
Index([u'1400.0', u'1500.0', u'1600.0', u'1700.0', u'1800.0', u'1900.0',
u'2000.0', u'2100.0', u'2200.0', u'2300.0', u'2400.0', u'2500.0',
u'2600.0', u'2700.0', u'2800.0', u'2900.0', u'3000.0', u'3100.0',
u'3200.0', u'3300.0', u'3400.0', u'3500.0', u'3600.0', u'3700.0',
u'3800.0', u'3900.0', u'4000.0', u'4100.0', u'4200.0', u'4300.0',
u'4400.0', u'4500.0', u'4600.0', u'4700.0', u'4800.0', u'4900.0',
u'5000.0', u'5100.0', u'5200.0', u'5300.0', u'5400.0', u'5500.0',
u'5600.0', u'5700.0', u'5800.0', u'5900.0', u'6000.0', u'6100.0',
u'6200.0', u'6300.0', u'6400.0', u'6500.0', u'6600.0', u'6700.0',
u'6800.0', u'6900.0', u'7000.0', u'7100.0', u'7200.0', u'7300.0',
u'7400.0', u'7500.0', u'7600.0', u'7700.0'],
dtype='object')
Index([u'1400.0', u'1500.0', u'1600.0', u'1700.0', u'1800.0', u'1900.0',
u'2000.0', u'2100.0', u'2200.0', u'2300.0', u'2400.0', u'2500.0',
u'2600.0', u'2700.0', u'2800.0', u'2900.0', u'3000.0', u'3100.0',
u'3200.0', u'3300.0', u'3400.0', u'3500.0', u'3600.0', u'3700.0',
u'3800.0', u'3900.0', u'4000.0', u'4100.0', u'4200.0', u'4300.0',
u'4400.0', u'4500.0', u'4600.0', u'4700.0', u'4800.0', u'4900.0',
u'5000.0', u'5100.0', u'5200.0', u'5300.0', u'5400.0', u'5500.0',
u'5600.0', u'5700.0', u'5800.0', u'5900.0', u'6000.0', u'6100.0',
u'6200.0', u'6300.0', u'6400.0', u'6500.0', u'6600.0', u'6700.0',
u'6800.0', u'6900.0', u'7000.0', u'7100.0', u'7200.0', u'7300.0',
u'7400.0', u'7500.0', u'7600.0', u'7700.0'],
dtype='object')
In [6]:
station_temp_hyps.data = station_temp_hyps.data - 273.15
station_temp_hyps.data.describe()
Out[6]:
1450
1550
1650
1750
1850
1950
2050
2150
2250
2350
...
6850
6950
7050
7150
7250
7350
7450
7550
7650
7750
count
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
...
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
mean
16.134795
15.484795
14.834795
14.184795
13.534795
12.884795
12.234795
11.584795
10.934795
10.284795
...
-18.918630
-19.568630
-20.218630
-20.868630
-21.518630
-22.168630
-22.818630
-23.468630
-24.118630
-24.768630
std
8.313342
8.313342
8.313342
8.313342
8.313342
8.313342
8.313342
8.313342
8.313342
8.313342
...
8.327385
8.327385
8.327385
8.327385
8.327385
8.327385
8.327385
8.327385
8.327385
8.327385
min
1.440000
0.790000
0.140000
-0.510000
-1.160000
-1.810000
-2.460000
-3.110000
-3.760000
-4.410000
...
-34.080000
-34.730000
-35.380000
-36.030000
-36.680000
-37.330000
-37.980000
-38.630000
-39.280000
-39.930000
25%
9.040000
8.390000
7.740000
7.090000
6.440000
5.790000
5.140000
4.490000
3.840000
3.190000
...
-26.280000
-26.930000
-27.580000
-28.230000
-28.880000
-29.530000
-30.180000
-30.830000
-31.480000
-32.130000
50%
16.440000
15.790000
15.140000
14.490000
13.840000
13.190000
12.540000
11.890000
11.240000
10.590000
...
-18.680000
-19.330000
-19.980000
-20.630000
-21.280000
-21.930000
-22.580000
-23.230000
-23.880000
-24.530000
75%
23.040000
22.390000
21.740000
21.090000
20.440000
19.790000
19.140000
18.490000
17.840000
17.190000
...
-12.180000
-12.830000
-13.480000
-14.130000
-14.780000
-15.430000
-16.080000
-16.730000
-17.380000
-18.030000
max
31.740000
31.090000
30.440000
29.790000
29.140000
28.490000
27.840000
27.190000
26.540000
25.890000
...
-2.880000
-3.530000
-4.180000
-4.830000
-5.480000
-6.130000
-6.780000
-7.430000
-8.080000
-8.730000
8 rows × 64 columns
In [7]:
corrected_temp_hyps.data.describe()
Out[7]:
1400.0
1500.0
1600.0
1700.0
1800.0
1900.0
2000.0
2100.0
2200.0
2300.0
...
6800.0
6900.0
7000.0
7100.0
7200.0
7300.0
7400.0
7500.0
7600.0
7700.0
count
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
...
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
mean
17.701452
17.297096
16.626904
16.030795
15.297233
14.784685
14.124438
13.697479
12.919836
12.103260
...
-18.306685
-19.333534
-19.494219
-20.070959
-20.651644
-20.629041
-21.348466
-22.113671
-21.373479
-24.797425
std
8.289469
8.335189
8.355992
8.393407
8.386968
8.382250
8.369800
8.358995
8.361076
8.355695
...
8.751661
8.764587
8.786613
8.784820
8.786821
8.785382
8.804826
8.819884
8.920145
8.799019
min
-0.940000
-1.350000
-2.010000
-2.640000
-3.310000
-3.790000
-4.380000
-4.770000
-5.490000
-6.240000
...
-35.410000
-36.610000
-36.790000
-37.450000
-38.140000
-38.070000
-38.910000
-39.800000
-38.990000
-42.860000
25%
11.000000
10.440000
9.580000
8.900000
8.140000
7.620000
6.990000
6.530000
5.680000
4.800000
...
-26.360000
-27.390000
-27.580000
-28.130000
-28.690000
-28.700000
-29.430000
-30.260000
-29.680000
-32.790000
50%
17.620000
17.170000
16.510000
15.920000
15.200000
14.760000
14.120000
13.760000
12.980000
12.150000
...
-18.230000
-19.250000
-19.410000
-20.010000
-20.590000
-20.590000
-21.280000
-22.040000
-21.360000
-24.650000
75%
25.230000
24.710000
24.260000
23.670000
22.880000
22.450000
21.840000
21.370000
20.560000
19.720000
...
-11.350000
-12.440000
-12.580000
-13.230000
-13.820000
-13.780000
-14.540000
-15.310000
-14.350000
-18.260000
max
35.380000
34.980000
34.280000
33.660000
32.870000
32.300000
31.580000
31.080000
30.260000
29.360000
...
-1.010000
-1.860000
-1.980000
-2.460000
-2.940000
-2.910000
-3.510000
-4.130000
-3.370000
-6.490000
8 rows × 64 columns
In [8]:
uncorrected_temp_hyps.data.describe()
Out[8]:
1400.0
1500.0
1600.0
1700.0
1800.0
1900.0
2000.0
2100.0
2200.0
2300.0
...
6800.0
6900.0
7000.0
7100.0
7200.0
7300.0
7400.0
7500.0
7600.0
7700.0
count
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.00000
365.000000
365.000000
365.000000
...
365.000000
365.000000
365.000000
365.000000
365.000000
365.000000
365.00000
365.000000
365.000000
365.000000
mean
17.000164
16.401068
15.623507
14.816438
14.063945
13.427205
12.76000
12.223315
11.430411
10.561644
...
-20.204055
-21.288740
-21.331890
-21.965562
-22.641534
-22.738767
-23.36211
-24.264932
-23.166466
-27.265699
std
7.522587
7.568178
7.589203
7.630471
7.626739
7.633017
7.62061
7.614093
7.614409
7.605798
...
7.947641
7.968682
7.976967
7.988620
7.999162
7.996180
8.01156
8.026520
8.045212
8.049260
min
-0.550000
-1.150000
-1.920000
-2.760000
-3.460000
-4.080000
-4.68000
-5.180000
-5.920000
-6.720000
...
-36.150000
-37.420000
-37.450000
-38.200000
-38.990000
-39.040000
-39.78000
-40.790000
-39.480000
-44.230000
25%
10.950000
10.250000
9.410000
8.520000
7.730000
7.050000
6.38000
5.820000
5.010000
4.130000
...
-27.770000
-28.770000
-28.850000
-29.440000
-30.160000
-30.250000
-30.93000
-31.790000
-30.740000
-34.620000
50%
17.320000
16.610000
15.810000
14.940000
14.160000
13.600000
12.96000
12.420000
11.580000
10.630000
...
-19.910000
-20.970000
-21.010000
-21.610000
-22.300000
-22.400000
-23.02000
-23.950000
-22.850000
-26.890000
75%
23.440000
22.950000
22.140000
21.470000
20.720000
20.260000
19.61000
19.180000
18.360000
17.450000
...
-13.790000
-14.950000
-15.000000
-15.690000
-16.440000
-16.550000
-17.23000
-18.250000
-17.000000
-21.390000
max
33.230000
32.650000
31.850000
31.030000
30.220000
29.540000
28.82000
28.220000
27.380000
26.440000
...
-4.380000
-5.270000
-5.290000
-5.810000
-6.370000
-6.460000
-6.97000
-7.730000
-6.740000
-10.320000
8 rows × 64 columns
In [9]:
station_temp_by_doy = (station_temp_hyps.data_by_doy() / len(station_temp_hyps.data.columns))
corrected_temp_by_doy = corrected_temp_hyps.data_by_doy() / len(corrected_temp_hyps.data.columns)
uncorrected_temp_by_doy = uncorrected_temp_hyps.data_by_doy() / len(uncorrected_temp_hyps.data.columns)
In [10]:
fig, ax = plt.subplots()
station_temp_by_doy.plot(color='k', linestyle='--', \
title=drainageID + " avg basin temperatures " + str(year), label='Stations')
corrected_temp_by_doy.plot(color='r', linestyle=':', label='Downscaled, Corrected')
uncorrected_temp_by_doy.plot(color='k', linestyle=':', label='Downscaled, Uncorrected')
ax.set_ylabel('Temperature')
legend = ax.legend()
fig.savefig("/Users/brodzik/charis/IN_Hunza_observed_vs_corrected_downscaled_temps.png")
In [ ]:
len(old_temp_hyps.data[str(3000 + 50)])
In [ ]:
elevations = [2000, 3000, 4000, 5000, 6000, 7000]
print old_temp_hyps.data
print temp_hyps.data
#print len(old_temp_hyps.data[str(elevations + 50)])
#print len(temp_hyps.data[str(float(elevations))][:-1])
In [ ]:
elevations = [2000, 3000, 4000, 5000, 6000, 7000]
nrows = 2
ncols = 3
fig, ax = plt.subplots(nrows, ncols)
fig.suptitle('IN_Hunza station vs. downscaled temperatures 2001', fontsize=18)
for i, el in enumerate(elevations):
row = i / ncols
col = i % ncols
this_ax = ax[row, col]
this_ax.scatter(old_temp_hyps.data[str(el + 50)], temp_hyps.data[str(float(el))])
lims = [
np.min([this_ax.get_xlim(), this_ax.get_ylim()]), # min of both axes
np.max([this_ax.get_xlim(), this_ax.get_ylim()]), # max of both axes
]
# now plot both limits against eachother
this_ax.plot(lims, lims, 'k-')
# this_ax.plot(lims, lims, 'k-', alpha=0.75, zorder=0)
this_ax.set_aspect('equal')
this_ax.set_xlim(lims)
this_ax.set_ylim(lims)
this_ax.set_xlabel('Stations')
this_ax.set_ylabel('Downscaled')
this_ax.set_title(str(el))
fig.tight_layout()
fig.savefig("IN_Hunza_observed_vs_downscaled_temps.scatter.png")
In [ ]:
station_temp_hyps.data.plot.box(vert=False,positions=np.arange(len(station_temp_hyps.data.columns)),
grid=True, title='Stations', label='stations')
uncorrected_temp_hyps.data.plot.box(vert=False, \
grid=True, title='Uncorrected dscale', \
color='k', \
positions=np.arange(len(uncorrected_temp_hyps.data.columns)),label='uncorrected_dscale')
corrected_temp_hyps.data.plot.box(vert=False, \
grid=True, title='Corrected dscale', \
color='r', \
positions=np.arange(len(corrected_temp_hyps.data.columns)),label='corrected_dscale')
#for i, el in enumerate(station_temp_hyps.data.columns):
# if i < 2:
# x = station_temp_hyps.data[el].values
# y.fill(int(el))
# ax.plot()
In [11]:
clim=(-10,10)
fig, ax = plt.subplots(1,5, figsize=(16,8))
ax[0].imshow(station_temp_hyps.data.values, cmap='RdBu_r', aspect='auto', clim=clim)
ax[0].set_title('Stations')
ax[1].imshow(corrected_temp_hyps.data.values, cmap='RdBu_r', aspect='auto', clim=clim)
ax[1].set_title('Dscale corrected')
ax[2].imshow(corrected_temp_hyps.data.values - station_temp_hyps.data.values, cmap='RdBu_r', aspect='auto', clim=clim)
ax[2].set_title('corrected-sta')
ax[3].imshow(uncorrected_temp_hyps.data.values, cmap='RdBu_r', aspect='auto', clim=clim)
ax[3].set_title('Dscale uncorrected')
ax[4].imshow(uncorrected_temp_hyps.data.values - station_temp_hyps.data.values, cmap='RdBu_r', aspect='auto', clim=clim)
ax[4].set_title('uncorrected-sta')
Out[11]:
<matplotlib.text.Text at 0x11773e710>
In [15]:
%cat /Users/brodzik/charis_training_2015/iPython_notebooks/temp_index_model/modis_tiles_config.ini
# Data processing configuration for MODIS-tile data
# used by CHARIS melt modeling
#
# Each section thereafter contains information for reading
# a type of data, including
#
# directory, filename pattern and expected data type (byte, float)
# Expected data type is only needed for flat binary data,
# since tif, h5 and nc file readers will take care of type issues
model_top_dir = ~/projects/CHARIS/
# Input data filename patterns
[input]
# Model forcing data that varies by time (annual, daily)
[[forcing]]
# Location of annual cubes of daily MODSCAG_GF (gap-filled modscag) snow data
[[[modscag_gf]]]
dir = %MODEL_TOP_DIR%/snow_cover/MODSCAG_GF/v05/%TILEID%/
pattern = MODSCAG_GF_Snow.v0.5.%TILEID%_%YYYY%.h5
# Location of annual cubes of daily grain size from
# MODSCAG_GF (gap-filled modscag) snow data
[[[grsize_scag]]]
dir = %MODEL_TOP_DIR%/grain_size/MODSCAG_GF/v05/%TILEID%/
pattern = MODSCAG_GF_GrainSize.v0.5.%TILEID%_%YYYY%.h5
# Location of downscaled temperature data
# "Corrected with station data:"
# dir = %MODEL_TOP_DIR%/forcing_data/Downscaled/ERA-Interim/day/corrected/v2/%TILEID%/%YYYY%/
# pattern = tsurf_%TILEID%_day_noinversion.%YYYYMMDD%.corrected.v2.bin
[[[temperature]]]
dir = %MODEL_TOP_DIR%/forcing_data/Downscaled/ERA-Interim/day/uncorrected/v1/%TILEID%/%YYYY%/
pattern = tsurf_%TILEID%_day_noinversion.%YYYYMMDD%.v1.bin
# Location of uncorrected downscaled ERA_interim temperature data
[[[uncorrected_temperature]]]
dir = %MODEL_TOP_DIR%/temperature/uncorrected/v0/%TILEID%/%YYYY%
pattern = tsurf_%TILEID%_day_%YYYYMMDD%.bin
# Location of 'corrected' downscaled ERA_interim temperature data
# 'corrected' temperatures are adjusted with station temperature
# climatologies at locations near available stations
[[[corrected_temperature]]]
dir = %MODEL_TOP_DIR%/temperature/corrected/v0/%TILEID%/%YYYY%
pattern = tsurf_%TILEID%_day_%YYYYMMDD%.corrected.v0.bin
[[fixed]]
[[[basin_mask]]]
dir = %MODEL_TOP_DIR%/basins/GDBD_subbasins/MODIS_tiles/
pattern = %DRAINAGEID%.basin_mask.%TILEID%.tif
[[[dem]]]
dir = %MODEL_TOP_DIR%/elevation_data/SRTMGL3/
pattern = SRTMGL3.v0.1.%TILEID%.tif
[[[modice_min05yr]]]
# Location of MODICE v04 5-year intersection data
dir = %MODEL_TOP_DIR%/snow_cover/modice.v0.4/min05yr_nc/
pattern = MODICE.v0.4.%TILEID%.%NSTRIKES%strike.min05yr.mask.nc
[hypsometry]
# MODICE by elevation
[[modice_min05yr_by_elevation]]
dir = %MODEL_TOP_DIR%/snow_cover/modice.v0.4/by_elevation
pattern = %DRAINAGEID%.%CONTOUR%.modicev04_%NSTRIKES%strike_area_by_elev.txt
[[snow_on_land_by_elevation]]
dir = %MODEL_TOP_DIR%/surface_area/by_elevation
pattern = %DRAINAGEID%.%YYYY%.%CONTOUR%.%SURFACE%_area_by_elev.txt
[[snow_on_ice_by_elevation]]
dir = %MODEL_TOP_DIR%/surface_area/by_elevation
pattern = %DRAINAGEID%.%YYYY%.%CONTOUR%.%ABLATION_METHOD%.%THRESHOLD%.%SURFACE%_area_by_elev.txt
[[exposed_glacier_ice_by_elevation]]
dir = %MODEL_TOP_DIR%/surface_area/by_elevation
pattern = %DRAINAGEID%.%YYYY%.%CONTOUR%.%ABLATION_METHOD%.%THRESHOLD%.%SURFACE%_area_by_elev.txt
# temperature by elevation
# dir = %MODEL_TOP_DIR%/forcing_data/Downscaled/ERA-Interim/day/corrected/by_elevation
# pattern = %DRAINAGEID%.%YYYY%.%CONTOUR%.corrected.v2.temperature_by_elev.txt
[[temperature_by_elevation]]
dir = %MODEL_TOP_DIR%/forcing_data/Downscaled/ERA-Interim/day/uncorrected/by_elevation
pattern = %DRAINAGEID%.%YYYY%.%CONTOUR%.v1.temperature_by_elev.txt
[[area_by_elevation]]
dir = %MODEL_TOP_DIR%/basins/GDBD_subbasins/by_elevation
pattern = %DRAINAGEID%.%CONTOUR%.area_by_elevation.txt
In [25]:
EGIFile = myEnv.hypsometry_filename(drainageID=drainageID, type='exposed_glacier_ice_by_elevation',year=2001,threshold=205)
EGIFile
Out[25]:
'/Users/brodzik/projects/CHARIS/surface_area/by_elevation/IN_Hunza_GDBD.2001.0100m.GRSIZE_SCAG.0205.exposed_glacier_ice_area_by_elev.txt'
In [19]:
SOLFile = myEnv.hypsometry_filename(drainageID=drainageID, type='snow_on_land_by_elevation', year=2001)
SOLFile
Out[19]:
'/Users/brodzik/projects/CHARIS/surface_area/by_elevation/IN_Hunza_GDBD.2001.0100m.snow_on_land_area_by_elev.txt'
In [21]:
SOIFile = myEnv.hypsometry_filename(drainageID=drainageID, type='snow_on_ice_by_elevation', year=2001,threshold=205)
SOIFile
Out[21]:
'/Users/brodzik/projects/CHARIS/surface_area/by_elevation/IN_Hunza_GDBD.2001.0100m.GRSIZE_SCAG.0205.snow_on_ice_area_by_elev.txt'
In [26]:
egi_hyps = Hypsometry(filename=EGIFile)
sol_hyps = Hypsometry(filename=SOLFile)
soi_hyps = Hypsometry(filename=SOIFile)
In [27]:
egi_hyps.data
Out[27]:
2900.0
3000.0
3100.0
3200.0
3300.0
3400.0
3500.0
3600.0
3700.0
3800.0
...
6800.0
6900.0
7000.0
7100.0
7200.0
7300.0
7400.0
7500.0
7600.0
7700.0
Date
2001-01-01
0.00
0.0
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.21
...
0.40
0.00
0.18
0.20
0.19
0.18
0.17
0.00
0.20
0.0
2001-01-02
0.00
0.0
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.21
...
0.40
0.00
0.00
0.18
0.19
0.18
0.17
0.00
0.18
0.0
2001-01-03
0.00
0.0
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
...
0.39
0.17
0.20
0.17
0.19
0.18
0.17
0.00
0.17
0.0
2001-01-04
0.00
0.0
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
...
0.38
0.18
0.20
0.00
0.19
0.18
0.18
0.00
0.17
0.0
2001-01-05
0.00
0.0
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
...
0.38
0.18
0.19
0.00
0.19
0.18
0.19
0.00
0.17
0.0
2001-01-06
0.00
0.0
0.00
0.00
0.35
0.00
0.00
0.00
0.00
0.00
...
0.38
0.17
0.19
0.00
0.18
0.19
0.20
0.00
0.00
0.0
2001-01-07
0.00
0.0
0.00
0.00
0.36
0.00
0.00
0.00
0.00
0.00
...
0.17
0.16
0.00
0.00
0.38
0.19
0.20
0.00
0.00
0.0
2001-01-08
0.00
0.0
0.00
0.00
0.35
0.00
0.00
0.00
0.00
0.00
...
0.16
0.15
0.00
0.19
0.38
0.19
0.37
0.00
0.00
0.0
2001-01-09
0.00
0.0
0.00
0.00
0.35
0.00
0.00
0.00
0.00
0.00
...
0.00
0.00
0.00
0.19
0.37
0.19
0.37
0.18
0.19
0.0
2001-01-10
0.00
0.0
0.00
0.00
0.34
0.15
0.00
0.00
0.00
0.00
...
0.00
0.00
0.00
0.55
0.55
0.18
0.19
0.18
0.20
0.0
2001-01-11
0.00
0.0
0.00
0.00
0.32
0.57
0.00
0.00
0.00
0.00
...
0.00
0.00
0.00
0.55
0.53
0.18
0.56
0.18
0.20
0.0
2001-01-12
0.00
0.0
0.15
0.00
0.32
0.44
0.16
0.00
0.00
0.00
...
0.17
0.00
0.00
0.38
0.51
0.17
0.37
0.18
0.19
0.0
2001-01-13
0.00
0.0
0.31
0.00
0.32
0.44
0.00
0.00
0.00
0.00
...
0.18
0.00
0.00
0.52
0.15
0.00
0.36
0.00
0.19
0.0
2001-01-14
0.14
0.0
0.48
0.18
0.32
0.44
0.00
0.00
0.00
0.00
...
0.19
0.00
0.18
0.70
0.13
0.00
0.68
0.14
0.19
0.0
2001-01-15
0.15
0.0
0.48
0.18
0.31
0.44
0.00
0.00
0.00
0.00
...
0.19
0.00
0.18
0.69
0.11
0.00
0.68
0.13
0.20
0.0
2001-01-16
0.15
0.0
0.47
0.18
0.29
0.27
0.00
0.00
0.00
0.00
...
0.19
0.16
0.18
0.55
0.00
0.00
0.67
0.14
0.20
0.0
2001-01-17
0.15
0.0
0.46
0.00
0.27
0.15
0.00
0.00
0.00
0.00
...
0.19
0.17
0.18
0.55
0.00
0.00
0.67
0.14
0.19
0.0
2001-01-18
0.14
0.0
0.45
0.00
0.00
0.00
0.00
0.00
0.15
0.16
...
0.37
0.36
0.18
0.55
0.00
0.00
0.17
0.14
0.19
0.0
2001-01-19
0.14
0.0
0.44
0.00
0.00
0.00
0.00
0.00
0.15
0.31
...
0.57
0.37
0.17
0.36
0.00
0.00
0.00
0.00
0.00
0.0
2001-01-20
0.13
0.0
0.42
0.00
0.00
0.00
0.00
0.00
0.14
0.30
...
0.57
0.37
0.17
0.36
0.00
0.00
0.00
0.00
0.00
0.0
2001-01-21
0.13
0.0
0.41
0.00
0.00
0.00
0.00
0.00
0.14
0.29
...
0.57
0.37
0.17
0.53
0.00
0.00
0.00
0.00
0.00
0.0
2001-01-22
0.00
0.0
0.41
0.00
0.00
0.00
0.00
0.00
0.14
0.60
...
0.72
0.55
0.00
0.53
0.13
0.00
0.00
0.00
0.00
0.0
2001-01-23
0.00
0.0
0.41
0.00
0.00
0.00
0.00
0.00
0.14
0.58
...
0.72
0.55
0.00
0.54
0.00
0.00
0.00
0.00
0.00
0.0
2001-01-24
0.00
0.0
0.42
0.15
0.00
0.00
0.00
0.00
0.14
0.57
...
0.55
0.55
0.00
0.54
0.00
0.00
0.00
0.00
0.00
0.0
2001-01-25
0.00
0.0
0.44
0.33
0.00
0.27
0.38
0.00
0.13
0.55
...
0.54
0.36
0.00
0.55
0.00
0.00
0.00
0.00
0.00
0.0
2001-01-26
0.00
0.0
0.49
0.50
0.24
0.41
0.47
0.00
0.13
0.66
...
0.52
0.18
0.00
0.55
0.00
0.00
0.00
0.00
0.00
0.0
2001-01-27
0.00
0.0
0.52
0.51
0.24
0.58
0.45
0.13
0.26
0.77
...
0.35
0.18
0.00
0.55
0.00
0.00
0.00
0.00
0.00
0.0
2001-01-28
0.00
0.0
0.52
0.88
0.66
0.55
0.43
0.13
0.52
0.86
...
0.36
0.18
0.00
0.37
0.00
0.00
0.00
0.00
0.00
0.0
2001-01-29
0.00
0.0
0.52
0.71
0.65
0.65
0.51
0.72
0.53
0.95
...
0.36
0.18
0.00
0.19
0.00
0.00
0.00
0.00
0.00
0.0
2001-01-30
0.00
0.0
0.65
0.56
0.64
0.63
0.52
1.00
0.55
1.17
...
0.36
0.19
0.00
0.19
0.00
0.00
0.00
0.00
0.00
0.0
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
2001-12-02
0.00
0.0
0.19
0.47
0.52
0.46
0.71
0.00
0.38
0.00
...
0.00
0.83
0.50
0.32
0.14
0.18
0.17
0.00
0.00
0.0
2001-12-03
0.00
0.0
0.16
0.09
0.29
0.38
0.58
0.00
0.40
0.00
...
0.36
0.67
0.30
0.13
0.14
0.00
0.17
0.00
0.00
0.0
2001-12-04
0.00
0.0
0.15
0.09
0.17
0.56
0.46
0.00
0.41
0.00
...
0.36
0.54
0.29
0.13
0.14
0.00
0.17
0.00
0.00
0.0
2001-12-05
0.00
0.0
0.17
0.10
0.16
0.56
0.46
0.00
0.41
0.00
...
0.50
0.55
0.28
0.13
0.00
0.00
0.00
0.00
0.00
0.0
2001-12-06
0.11
0.0
0.32
0.12
0.14
0.41
0.15
0.00
0.40
0.00
...
0.82
0.39
0.43
0.12
0.13
0.00
0.00
0.00
0.00
0.0
2001-12-07
0.12
0.0
0.38
0.00
0.12
0.36
0.13
0.00
0.20
0.00
...
0.93
0.58
0.62
0.12
0.15
0.00
0.00
0.00
0.00
0.0
2001-12-08
0.12
0.0
0.13
0.00
0.11
0.22
0.00
0.00
0.20
0.21
...
1.07
0.60
0.77
0.32
0.16
0.00
0.00
0.00
0.00
0.0
2001-12-09
0.12
0.0
0.14
0.00
0.11
0.12
0.21
0.00
0.20
0.20
...
1.01
0.61
0.61
0.50
0.56
0.00
0.00
0.00
0.00
0.0
2001-12-10
0.13
0.0
0.15
0.00
0.14
0.12
0.43
0.00
0.20
0.20
...
1.24
0.61
0.44
0.89
0.58
0.00
0.19
0.00
0.20
0.0
2001-12-11
0.13
0.0
0.15
0.00
0.17
0.14
1.06
0.00
0.20
0.20
...
1.49
0.59
0.64
1.09
0.80
0.18
0.19
0.00
0.20
0.0
2001-12-12
0.14
0.0
0.15
0.00
0.00
0.00
1.03
0.00
0.19
0.20
...
1.48
0.75
0.66
1.10
0.81
0.19
0.39
0.00
0.19
0.0
2001-12-13
0.00
0.0
0.00
0.00
0.00
0.00
0.97
0.00
0.33
0.19
...
1.43
0.37
0.72
0.75
0.42
0.00
0.20
0.00
0.17
0.0
2001-12-14
0.00
0.0
0.00
0.00
0.00
0.00
0.89
0.00
0.29
0.18
...
1.58
0.39
0.87
0.55
0.41
0.18
0.21
0.00
0.15
0.0
2001-12-15
0.00
0.0
0.00
0.00
0.00
0.00
0.85
0.00
0.27
0.18
...
1.29
0.39
0.89
0.68
0.41
0.18
0.21
0.00
0.00
0.0
2001-12-16
0.00
0.0
0.00
0.00
0.00
0.00
0.82
0.14
0.27
0.17
...
1.47
0.19
0.71
0.49
0.41
0.17
0.21
0.16
0.00
0.0
2001-12-17
0.00
0.0
0.00
0.00
0.00
0.00
0.33
0.00
0.14
0.16
...
1.47
0.59
0.70
0.47
0.21
0.16
0.21
0.00
0.00
0.0
2001-12-18
0.00
0.0
0.00
0.00
0.00
0.00
0.33
0.00
0.15
0.15
...
1.16
0.59
0.70
0.31
0.21
0.16
0.21
0.00
0.00
0.0
2001-12-19
0.00
0.0
0.00
0.00
0.00
0.00
0.34
0.00
0.00
0.00
...
0.96
0.77
0.88
0.30
0.41
0.15
0.20
0.20
0.00
0.0
2001-12-20
0.00
0.0
0.00
0.00
0.00
0.00
0.34
0.00
0.00
0.00
...
1.26
0.78
1.03
0.29
0.21
0.15
0.00
0.00
0.00
0.0
2001-12-21
0.00
0.0
0.00
0.00
0.00
0.00
0.34
0.00
0.00
0.00
...
1.25
0.95
0.85
0.46
0.20
0.16
0.00
0.00
0.00
0.0
2001-12-22
0.00
0.0
0.00
0.00
0.00
0.00
0.34
0.00
0.00
0.00
...
1.03
0.93
0.86
0.46
0.18
0.00
0.00
0.00
0.00
0.0
2001-12-23
0.00
0.0
0.00
0.00
0.00
0.00
0.34
0.00
0.00
0.00
...
1.01
1.10
1.25
0.46
0.37
0.38
0.19
0.00
0.00
0.0
2001-12-24
0.00
0.0
0.00
0.00
0.00
0.00
0.34
0.00
0.00
0.00
...
1.36
1.07
1.37
0.46
0.71
0.37
0.19
0.00
0.00
0.0
2001-12-25
0.00
0.0
0.00
0.00
0.00
0.00
0.33
0.00
0.00
0.00
...
1.56
0.89
1.38
0.46
0.70
0.37
0.19
0.00
0.00
0.0
2001-12-26
0.00
0.0
0.00
0.00
0.00
0.00
0.33
0.00
0.00
0.00
...
1.56
0.88
1.38
0.64
0.69
0.36
0.20
0.00
0.00
0.0
2001-12-27
0.00
0.0
0.00
0.00
0.17
0.00
0.32
0.00
0.00
0.00
...
1.27
0.70
0.82
0.48
0.71
0.36
0.20
0.00
0.18
0.0
2001-12-28
0.00
0.0
0.00
0.00
0.17
0.16
0.70
0.00
0.31
0.30
...
1.31
0.73
0.73
0.91
0.93
0.19
0.21
0.00
0.19
0.0
2001-12-29
0.00
0.0
0.00
0.00
0.17
0.32
0.62
0.12
0.29
0.32
...
1.29
1.13
0.73
0.74
1.48
0.18
0.20
0.00
0.00
0.0
2001-12-30
0.00
0.0
0.00
0.00
0.15
0.39
0.66
0.12
0.25
0.32
...
1.13
1.28
0.71
0.54
0.91
0.17
0.20
0.00
0.00
0.0
2001-12-31
0.00
0.0
0.00
0.00
0.29
0.29
0.70
0.12
0.14
0.32
...
0.69
1.01
0.84
0.49
0.86
0.15
0.39
0.00
0.20
0.0
365 rows × 49 columns
In [28]:
sol_hyps.data
Out[28]:
1400.0
1500.0
1600.0
1700.0
1800.0
1900.0
2000.0
2100.0
2200.0
2300.0
...
5800.0
5900.0
6000.0
6100.0
6200.0
6300.0
6400.0
6500.0
6600.0
6700.0
Date
2001-01-01
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.09
0.22
1.15
...
1.15
0.20
0.64
0.02
0.37
0.52
0.44
0.22
0.05
0.06
2001-01-02
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.09
0.22
1.15
...
1.14
0.21
0.64
0.02
0.38
0.49
0.45
0.21
0.05
0.06
2001-01-03
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.13
0.21
1.18
...
1.14
0.21
0.65
0.02
0.39
0.46
0.46
0.21
0.05
0.07
2001-01-04
0.0
0.0
0.0
0.0
0.0
0.02
0.03
0.14
0.21
1.16
...
1.14
0.23
0.65
0.02
0.41
0.45
0.50
0.21
0.03
0.08
2001-01-05
0.0
0.0
0.0
0.0
0.0
0.02
0.04
0.15
0.21
1.14
...
1.14
0.25
0.66
0.02
0.43
0.45
0.54
0.21
0.03
0.08
2001-01-06
0.0
0.0
0.0
0.0
0.0
0.02
0.00
0.17
0.20
1.11
...
1.13
0.27
0.65
0.02
0.46
0.45
0.60
0.22
0.03
0.08
2001-01-07
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.15
0.20
1.03
...
1.12
0.32
0.66
0.02
0.47
0.45
0.65
0.23
0.05
0.09
2001-01-08
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.14
0.20
0.90
...
1.12
0.34
0.67
0.02
0.45
0.46
0.69
0.23
0.05
0.09
2001-01-09
0.0
0.0
0.0
0.0
0.0
0.00
0.04
0.09
0.17
0.80
...
1.12
0.35
0.70
0.03
0.41
0.48
0.71
0.23
0.06
0.09
2001-01-10
0.0
0.0
0.0
0.0
0.0
0.00
0.04
0.06
0.16
0.81
...
1.11
0.34
0.73
0.03
0.39
0.49
0.73
0.23
0.07
0.09
2001-01-11
0.0
0.0
0.0
0.0
0.0
0.00
0.04
0.06
0.14
0.79
...
1.09
0.35
0.77
0.03
0.39
0.49
0.74
0.23
0.07
0.09
2001-01-12
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.02
0.09
0.73
...
1.09
0.34
0.79
0.03
0.38
0.49
0.73
0.22
0.07
0.09
2001-01-13
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.00
0.07
0.57
...
1.10
0.35
0.78
0.03
0.37
0.48
0.67
0.22
0.08
0.10
2001-01-14
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.00
0.04
0.53
...
1.13
0.35
0.77
0.03
0.37
0.47
0.64
0.22
0.09
0.11
2001-01-15
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.00
0.04
0.50
...
1.15
0.35
0.78
0.03
0.38
0.48
0.62
0.22
0.09
0.12
2001-01-16
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.00
0.04
0.47
...
1.18
0.36
0.79
0.03
0.39
0.49
0.63
0.22
0.10
0.12
2001-01-17
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.04
0.07
0.57
...
1.20
0.37
0.80
0.04
0.41
0.50
0.63
0.22
0.10
0.12
2001-01-18
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.06
0.08
0.62
...
1.22
0.38
0.80
0.04
0.43
0.52
0.64
0.22
0.10
0.12
2001-01-19
0.0
0.0
0.0
0.0
0.0
0.00
0.03
0.08
0.11
0.67
...
1.25
0.38
0.81
0.04
0.46
0.53
0.65
0.22
0.10
0.12
2001-01-20
0.0
0.0
0.0
0.0
0.0
0.00
0.05
0.08
0.11
0.71
...
1.28
0.39
0.81
0.04
0.48
0.56
0.66
0.22
0.10
0.12
2001-01-21
0.0
0.0
0.0
0.0
0.0
0.00
0.04
0.08
0.09
0.71
...
1.30
0.40
0.82
0.05
0.49
0.58
0.66
0.22
0.10
0.12
2001-01-22
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.04
0.08
0.60
...
1.32
0.40
0.82
0.05
0.51
0.59
0.67
0.22
0.09
0.12
2001-01-23
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.02
0.03
0.53
...
1.34
0.40
0.82
0.05
0.52
0.61
0.69
0.22
0.09
0.12
2001-01-24
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.02
0.03
0.50
...
1.36
0.41
0.83
0.07
0.52
0.61
0.70
0.22
0.09
0.12
2001-01-25
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.02
0.03
0.53
...
1.36
0.41
0.85
0.07
0.52
0.61
0.70
0.22
0.10
0.12
2001-01-26
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.02
0.03
0.49
...
1.33
0.40
0.85
0.07
0.51
0.59
0.71
0.21
0.10
0.13
2001-01-27
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.02
0.03
0.36
...
1.28
0.38
0.87
0.07
0.50
0.57
0.72
0.21
0.10
0.13
2001-01-28
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.05
0.03
0.20
...
1.22
0.36
0.87
0.07
0.48
0.56
0.73
0.21
0.10
0.13
2001-01-29
0.0
0.0
0.0
0.0
0.0
0.02
0.02
0.03
0.02
0.21
...
1.18
0.35
0.88
0.06
0.50
0.56
0.75
0.21
0.10
0.13
2001-01-30
0.0
0.0
0.0
0.0
0.0
0.02
0.00
0.03
0.00
0.17
...
1.17
0.35
0.89
0.06
0.50
0.55
0.77
0.21
0.10
0.13
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
2001-12-02
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.06
0.00
0.40
...
1.68
0.59
1.43
0.32
0.56
0.65
0.83
0.25
0.26
0.26
2001-12-03
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.06
0.00
0.33
...
1.67
0.55
1.40
0.24
0.56
0.64
0.82
0.25
0.25
0.25
2001-12-04
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.06
0.00
0.28
...
1.65
0.52
1.37
0.18
0.57
0.64
0.81
0.25
0.25
0.24
2001-12-05
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.06
0.00
0.35
...
1.63
0.52
1.34
0.15
0.58
0.65
0.82
0.25
0.25
0.25
2001-12-06
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.05
0.00
0.47
...
1.60
0.55
1.31
0.16
0.59
0.68
0.84
0.26
0.26
0.27
2001-12-07
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.03
0.00
0.53
...
1.58
0.55
1.27
0.17
0.59
0.72
0.85
0.26
0.27
0.29
2001-12-08
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.03
0.00
0.62
...
1.57
0.54
1.23
0.19
0.60
0.74
0.85
0.27
0.28
0.30
2001-12-09
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.03
0.00
0.68
...
1.55
0.50
1.17
0.22
0.59
0.75
0.82
0.28
0.27
0.30
2001-12-10
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.02
0.00
0.68
...
1.52
0.43
1.09
0.26
0.57
0.74
0.75
0.27
0.25
0.29
2001-12-11
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.00
0.00
0.71
...
1.49
0.37
1.02
0.29
0.55
0.72
0.68
0.27
0.22
0.28
2001-12-12
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.04
0.00
0.78
...
1.47
0.32
0.98
0.32
0.53
0.70
0.63
0.27
0.20
0.27
2001-12-13
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.04
0.03
1.07
...
1.41
0.30
0.93
0.35
0.53
0.68
0.60
0.27
0.19
0.27
2001-12-14
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.06
0.06
1.35
...
1.35
0.26
0.89
0.37
0.54
0.66
0.60
0.27
0.19
0.26
2001-12-15
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.08
0.10
1.47
...
1.34
0.26
0.88
0.39
0.56
0.67
0.61
0.27
0.19
0.26
2001-12-16
0.0
0.0
0.0
0.0
0.0
0.00
0.04
0.08
0.13
1.49
...
1.35
0.25
0.87
0.39
0.59
0.69
0.63
0.28
0.20
0.26
2001-12-17
0.0
0.0
0.0
0.0
0.0
0.00
0.04
0.06
0.09
1.46
...
1.39
0.26
0.88
0.39
0.62
0.71
0.66
0.28
0.21
0.26
2001-12-18
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.05
0.09
1.37
...
1.40
0.26
0.89
0.40
0.66
0.74
0.70
0.28
0.23
0.27
2001-12-19
0.0
0.0
0.0
0.0
0.0
0.00
0.02
0.02
0.11
1.33
...
1.44
0.26
0.91
0.40
0.72
0.78
0.74
0.30
0.26
0.28
2001-12-20
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.00
0.08
1.23
...
1.46
0.27
0.94
0.40
0.80
0.82
0.79
0.33
0.30
0.30
2001-12-21
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.00
0.06
1.16
...
1.44
0.28
0.96
0.40
0.83
0.84
0.82
0.34
0.32
0.31
2001-12-22
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.00
0.04
1.16
...
1.42
0.29
1.00
0.41
0.83
0.84
0.82
0.33
0.31
0.30
2001-12-23
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.02
0.04
1.12
...
1.37
0.31
1.03
0.41
0.81
0.84
0.81
0.31
0.29
0.29
2001-12-24
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.02
0.06
1.10
...
1.31
0.32
1.06
0.41
0.78
0.83
0.79
0.28
0.27
0.27
2001-12-25
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.03
0.04
1.06
...
1.26
0.33
1.09
0.41
0.76
0.82
0.77
0.27
0.26
0.26
2001-12-26
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.03
0.04
0.97
...
1.20
0.35
1.11
0.39
0.73
0.80
0.74
0.25
0.24
0.24
2001-12-27
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.02
0.02
0.90
...
1.16
0.37
1.11
0.37
0.71
0.77
0.72
0.25
0.24
0.23
2001-12-28
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.02
0.02
0.77
...
1.16
0.41
1.12
0.33
0.68
0.72
0.70
0.25
0.25
0.24
2001-12-29
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.00
0.02
0.61
...
1.22
0.46
1.18
0.34
0.67
0.68
0.70
0.26
0.26
0.27
2001-12-30
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.00
0.00
0.56
...
1.28
0.50
1.25
0.37
0.67
0.66
0.71
0.28
0.29
0.31
2001-12-31
0.0
0.0
0.0
0.0
0.0
0.00
0.00
0.00
0.00
0.47
...
1.29
0.52
1.31
0.36
0.63
0.64
0.74
0.29
0.30
0.32
365 rows × 54 columns
In [29]:
soi_hyps.data
Out[29]:
2900.0
3000.0
3100.0
3200.0
3300.0
3400.0
3500.0
3600.0
3700.0
3800.0
...
6800.0
6900.0
7000.0
7100.0
7200.0
7300.0
7400.0
7500.0
7600.0
7700.0
Date
2001-01-01
0.18
0.38
0.94
2.26
2.82
3.57
4.40
5.63
4.71
6.71
...
13.31
11.80
7.74
8.20
6.18
3.14
2.26
2.44
0.84
0.15
2001-01-02
0.17
0.38
0.88
2.21
2.72
3.54
4.30
5.43
4.65
6.48
...
13.29
11.66
7.91
8.26
6.15
3.15
2.20
2.37
0.83
0.14
2001-01-03
0.15
0.36
0.84
2.16
2.66
3.53
4.23
5.26
4.58
6.53
...
13.28
11.45
7.73
8.37
6.15
3.15
2.17
2.29
0.83
0.13
2001-01-04
0.14
0.33
0.80
2.14
2.68
3.57
4.27
5.21
4.54
6.52
...
13.32
11.48
7.76
8.58
6.19
3.17
2.16
2.28
0.83
0.13
2001-01-05
0.14
0.32
0.79
2.12
2.73
3.63
4.35
5.33
4.55
6.62
...
13.39
11.58
7.80
8.59
6.21
3.19
2.19
2.31
0.84
0.13
2001-01-06
0.14
0.31
0.81
2.13
2.43
3.68
4.45
5.56
4.63
6.81
...
13.46
11.69
7.86
8.58
6.25
3.20
2.21
2.36
1.02
0.14
2001-01-07
0.14
0.32
0.85
2.17
2.47
3.74
4.51
5.76
4.72
6.99
...
13.68
11.79
8.07
8.55
6.07
3.21
2.24
2.42
1.03
0.15
2001-01-08
0.15
0.32
0.88
2.25
2.51
3.76
4.51
5.84
4.78
7.07
...
13.65
11.84
8.06
8.28
6.06
3.20
2.09
2.48
1.04
0.15
2001-01-09
0.15
0.33
0.89
2.29
2.51
3.77
4.46
5.80
4.78
7.10
...
13.69
11.97
8.00
8.17
6.06
3.17
2.10
2.35
0.85
0.16
2001-01-10
0.14
0.34
0.87
2.28
2.50
3.56
4.40
5.71
4.73
7.09
...
13.52
11.90
7.95
7.70
5.86
3.13
2.28
2.37
0.83
0.16
2001-01-11
0.12
0.35
0.84
2.24
2.48
3.06
4.34
5.59
4.63
7.06
...
13.37
11.84
7.91
7.62
5.80
3.10
1.89
2.35
0.82
0.17
2001-01-12
0.11
0.35
0.68
2.20
2.46
3.21
4.21
5.58
4.62
7.09
...
13.14
11.83
7.94
7.76
5.73
3.10
2.05
2.32
0.80
0.17
2001-01-13
0.12
0.34
0.52
2.17
2.42
3.26
4.48
5.66
4.70
7.16
...
13.27
11.84
8.00
7.59
6.00
3.31
2.05
2.49
0.80
0.17
2001-01-14
0.00
0.33
0.34
1.98
2.37
3.25
4.60
5.71
4.78
7.19
...
13.49
11.88
7.89
7.40
5.96
3.38
1.73
2.37
0.80
0.18
2001-01-15
0.00
0.33
0.34
1.97
2.36
3.23
4.66
5.71
4.82
7.17
...
13.62
11.89
7.92
7.41
5.96
3.42
1.74
2.38
0.80
0.18
2001-01-16
0.00
0.33
0.34
1.95
2.37
3.39
4.66
5.65
4.79
7.13
...
13.65
11.74
7.92
7.59
6.10
3.43
1.75
2.37
0.79
0.18
2001-01-17
0.00
0.33
0.34
2.08
2.39
3.52
4.65
5.60
4.75
7.09
...
13.67
11.74
7.93
7.64
6.12
3.43
1.75
2.38
0.79
0.18
2001-01-18
0.00
0.32
0.33
2.03
2.65
3.68
4.62
5.57
4.58
6.88
...
13.48
11.54
7.91
7.68
6.16
3.43
2.23
2.39
0.79
0.18
2001-01-19
0.00
0.32
0.33
1.99
2.66
3.68
4.59
5.55
4.58
6.70
...
13.26
11.54
7.91
7.94
6.20
3.43
2.40
2.53
0.97
0.18
2001-01-20
0.00
0.31
0.33
1.95
2.66
3.69
4.57
5.54
4.59
6.69
...
13.24
11.54
7.90
8.00
6.24
3.43
2.39
2.54
0.97
0.18
2001-01-21
0.00
0.31
0.32
1.91
2.65
3.70
4.55
5.53
4.59
6.67
...
13.23
11.53
7.88
7.90
6.27
3.43
2.38
2.55
0.96
0.18
2001-01-22
0.12
0.30
0.32
1.90
2.66
3.72
4.56
5.53
4.59
6.32
...
13.08
11.35
8.05
7.94
6.17
3.43
2.38
2.55
0.96
0.18
2001-01-23
0.12
0.30
0.32
1.92
2.68
3.73
4.56
5.48
4.54
6.28
...
13.07
11.34
8.04
7.95
6.34
3.41
2.37
2.55
0.95
0.17
2001-01-24
0.11
0.30
0.32
1.76
2.61
3.69
4.48
5.33
4.45
6.14
...
13.21
11.29
8.00
7.93
6.33
3.41
2.37
2.54
0.94
0.17
2001-01-25
0.10
0.29
0.32
1.51
2.43
3.27
3.91
5.03
4.30
5.88
...
13.15
11.45
7.96
7.90
6.32
3.41
2.35
2.54
0.94
0.17
2001-01-26
0.10
0.28
0.31
1.31
2.10
3.02
3.70
4.77
4.21
5.51
...
13.09
11.60
7.93
7.87
6.31
3.41
2.34
2.54
0.95
0.17
2001-01-27
0.10
0.26
0.30
1.33
2.10
2.80
3.69
4.52
4.05
5.26
...
13.21
11.58
7.92
7.86
6.33
3.41
2.33
2.54
0.96
0.16
2001-01-28
0.10
0.25
0.29
0.95
1.64
2.75
3.72
4.49
3.78
5.09
...
13.18
11.58
7.91
8.05
6.33
3.41
2.33
2.54
0.97
0.16
2001-01-29
0.10
0.24
0.28
1.08
1.59
2.55
3.60
3.95
3.75
4.93
...
13.21
11.63
7.93
8.23
6.36
3.42
2.32
2.54
0.98
0.16
2001-01-30
0.10
0.23
0.14
1.20
1.59
2.50
3.61
3.78
3.72
4.75
...
13.27
11.70
7.98
8.25
6.39
3.42
2.33
2.55
0.98
0.16
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
2001-12-02
0.09
0.21
0.33
1.02
1.76
2.75
3.42
4.87
3.79
6.18
...
14.33
10.97
7.91
8.05
6.10
3.31
2.17
2.45
0.99
0.17
2001-12-03
0.11
0.23
0.33
1.46
2.10
2.92
3.74
5.11
3.95
6.46
...
13.99
11.17
8.08
8.27
6.09
3.50
2.16
2.45
0.98
0.17
2001-12-04
0.11
0.24
0.34
1.52
2.34
2.84
4.04
5.34
4.12
6.69
...
13.97
11.33
8.06
8.29
6.08
3.49
2.16
2.45
0.97
0.17
2001-12-05
0.12
0.26
0.37
1.55
2.36
2.90
4.14
5.49
4.24
6.82
...
13.76
11.33
8.06
8.32
6.24
3.48
2.33
2.45
0.97
0.17
2001-12-06
0.00
0.29
0.30
1.64
2.39
3.09
4.49
5.56
4.31
6.87
...
13.39
11.53
7.90
8.34
6.21
3.47
2.34
2.46
0.97
0.17
2001-12-07
0.00
0.33
0.32
1.97
2.47
3.15
4.49
5.56
4.52
6.90
...
13.36
11.44
7.73
8.40
6.35
3.48
2.37
2.48
1.00
0.17
2001-12-08
0.00
0.34
0.61
2.04
2.50
3.26
4.59
5.60
4.52
6.70
...
13.38
11.53
7.64
8.25
6.41
3.49
2.38
2.50
1.03
0.17
2001-12-09
0.00
0.32
0.60
2.00
2.48
3.37
4.38
5.68
4.53
6.71
...
13.64
11.61
7.86
8.14
6.07
3.51
2.38
2.51
1.05
0.17
2001-12-10
0.00
0.30
0.59
1.95
2.48
3.43
4.26
5.76
4.56
6.74
...
13.62
11.70
8.07
7.79
6.10
3.52
2.19
2.51
0.87
0.17
2001-12-11
0.00
0.28
0.57
1.91
2.48
3.53
3.74
5.82
4.58
6.77
...
13.53
11.77
7.91
7.61
5.91
3.35
2.19
2.51
0.88
0.18
2001-12-12
0.00
0.28
0.57
1.90
2.63
3.69
3.76
5.79
4.53
6.74
...
13.54
11.67
7.93
7.59
5.91
3.34
2.02
2.50
0.87
0.18
2001-12-13
0.14
0.29
0.75
1.92
2.56
3.62
3.68
5.64
4.24
6.60
...
13.41
12.08
7.92
7.89
6.26
3.51
2.23
2.49
0.85
0.18
2001-12-14
0.14
0.31
0.78
1.94
2.48
3.53
3.58
5.47
4.11
6.43
...
13.04
12.05
7.80
8.03
6.22
3.30
2.26
2.46
0.83
0.17
2001-12-15
0.15
0.32
0.80
1.93
2.41
3.46
3.51
5.36
4.04
6.30
...
13.27
12.06
7.79
7.88
6.19
3.27
2.28
2.45
0.95
0.17
2001-12-16
0.15
0.33
0.83
1.95
2.41
3.41
3.51
5.21
4.04
6.27
...
13.11
12.30
7.97
8.08
6.17
3.26
2.31
2.29
0.94
0.17
2001-12-17
0.15
0.35
0.86
2.02
2.53
3.38
4.01
5.46
4.23
6.37
...
13.23
12.04
7.98
8.10
6.38
3.26
2.35
2.48
0.93
0.17
2001-12-18
0.15
0.36
0.87
2.09
2.68
3.36
4.02
5.61
4.35
6.56
...
13.66
12.15
7.97
8.29
6.38
3.27
2.37
2.49
0.93
0.17
2001-12-19
0.15
0.36
0.89
2.09
2.73
3.31
4.03
5.58
4.64
6.80
...
13.94
11.95
7.75
8.27
6.16
3.30
2.37
2.27
0.91
0.17
2001-12-20
0.15
0.34
0.91
2.05
2.70
3.27
4.04
5.39
4.82
6.79
...
13.66
11.83
7.50
8.23
6.32
3.32
2.54
2.46
0.90
0.17
2001-12-21
0.15
0.33
0.93
2.05
2.67
3.26
4.05
5.28
4.89
6.77
...
13.65
11.55
7.60
8.03
6.32
3.33
2.52
2.45
0.90
0.17
2001-12-22
0.16
0.32
0.93
2.05
2.66
3.25
4.04
5.23
4.86
6.77
...
13.86
11.53
7.58
8.04
6.35
3.50
2.53
2.45
0.92
0.17
2001-12-23
0.16
0.33
0.94
2.05
2.64
3.26
4.05
5.21
4.81
6.77
...
13.89
11.33
7.23
8.07
6.20
3.14
2.34
2.45
0.93
0.17
2001-12-24
0.17
0.33
0.94
2.04
2.63
3.28
4.06
5.21
4.75
6.77
...
13.54
11.36
7.14
8.11
5.91
3.17
2.36
2.46
0.94
0.18
2001-12-25
0.17
0.34
0.94
2.05
2.64
3.31
4.08
5.21
4.69
6.81
...
13.36
11.57
7.17
8.14
5.99
3.20
2.37
2.50
0.96
0.18
2001-12-26
0.18
0.35
0.94
2.05
2.67
3.35
4.10
5.25
4.64
6.87
...
13.34
11.62
7.20
8.03
6.06
3.23
2.39
2.56
0.97
0.18
2001-12-27
0.18
0.35
0.94
2.05
2.56
3.41
4.15
5.29
4.62
6.91
...
13.61
11.88
7.76
8.24
6.11
3.27
2.40
2.61
0.80
0.18
2001-12-28
0.17
0.36
0.94
2.04
2.58
3.32
3.76
5.27
4.29
6.61
...
13.54
11.92
7.86
7.88
5.93
3.46
2.38
2.63
0.78
0.18
2001-12-29
0.17
0.35
0.92
2.00
2.55
3.23
3.70
5.15
4.30
6.56
...
13.55
11.49
7.84
8.09
5.34
3.47
2.32
2.59
0.95
0.17
2001-12-30
0.16
0.35
0.90
1.97
2.50
3.18
3.45
5.16
4.30
6.57
...
13.72
11.26
7.83
8.28
5.84
3.46
2.33
2.54
0.94
0.16
2001-12-31
0.16
0.35
0.90
1.97
2.33
3.26
3.31
5.24
4.38
6.60
...
14.11
11.39
7.66
8.21
5.79
3.40
2.16
2.48
0.76
0.17
365 rows × 49 columns
In [ ]:
Content source: mjbrodzik/ipython_notebooks
Similar notebooks: